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Motivated by recent discussion on possible quantum critical behavior in the coupled Cu- 
tetrahedra system Cu2Te2 0sBr2, we present a comparative ab initio study of the electronic prop- 
erties of Cu2Te2C>5Br2 and the isostructural Cu2Te205Ci2. A detailed investigation of the copper- 
copper interaction pathes reveals that the halogen-ions play an important role in the inter-tetrahedral 
couplings via X4-rings (X=Br, CI). We find that, contrary to initial indications, both systems show 
a similar electronic behavior with long range exchange pathes mediated by the X4-rings. 
PACS numbers: 75.30.Gw, 75.10.Jm, 78.30.-j 



The recently discovered [1] spin-tetrahedral com- 
pounds Cu 2 Te 2 05X 2 (X=C1, Br) open the possibility to 
study the interplay between localized many-body tetra- 
hedral cluster excitations and inter-tetrahedral magnetic 
couplings leading to a quantum-phase transition with 
various possible ordered states [2] . Transitions to ordered 
states have been observed experimentally [3] in these 
compounds with T^ Br) = 11.4 K and T^ C1) = 18.2 K. 
These phase-transitions exhibit unusual magnetic-field 
dependences which have been linked to the closeness to 
a quantum-critical point [3,4]. Unconventional Raman- 
scattering has been found in the magnetic channel [2,3] 
and the occurrence of low-lying singlet excitations has 
been proposed [1] and observed by Raman together with 
a longitudinal magnon [4] . 

The nature of the ordered states in Cu2Te 2 05X2 has 
not yet been definitely settled. There is some evidence 
for a Neel state from thermodynamic and susceptibility 
experiments and Cu2Te205Br2 has been proposed to be 
closer to a non-magnetic singlet state than Cu2Te20sCl2 
[4]. Thus, it is important to examine the microscopic 
behavior of both Cu2Te20sBr2 and Cu2Te20sCl2 by ab- 
initio methods and investigate whether the electronic 
properties reveal some non-trivial differences between the 
two systems. The magnetic exchange coupling parame- 
ters in these systems estimated from susceptibility mea- 
surements [1] are small, of the order 40 — 50 K, implying 
a small bandwith for the electronically active Cu-3g? or- 
bitals close to the Fermi-edge. 

Here we present a comprehensive first-principles den- 
sity functional theory (DFT) study of the electronic 
properties of these systems within the local spin den- 
sity (LSDA) and the generalized gradient approximation 
(GGA) [5]. 

Crystal structure.- Both Cu2Te20sBr2 and 
Cu2Te2 0sCl2 systems crystallize in the non- 
centrosymmetric tetragonal PA space group with two 
formula units per unit cell. Cu2Te205Br2 with lattice 
parameters a=7.84 A, c=6.38 A has a larger unit cell 



than Cu2Te 2 5 Cl2 with a=7.62 A, c=6.32 A. 

In Fig. 1 we show the crystal structure for 
Cu2Te2C>5X2. The four equivalent Cu 2+ -ions in the unit 
cell of these systems form distorted tetrahedra which 
are located on a tetragonal lattice. The structure con- 
tains three inequivalent oxygen positions: The intra- 
tetrahedral 01, and the inter-tetrahedral 02 and 03. 
Of special interest are the four equivalent halogen-ions, 
X, per unit-cell which form together with the 01 and 
02 the X-Ol-Ol-02 distorted square with the Cu 2+ - 
ions in their respective centers, see inset of Fig. 1. 
The inter-tetrahedron distances are slightly smaller in 
Cu2Te20sCl2 than in Cu2Te205Br2 due to the size dif- 
ference of the respective halogen-ions. The Cu-Cu intra- 
tetrahedron distances are, on the other hand, slightly 
larger in the Cl-compound. 

Band-structure.- In Fig. 2 we present the non- 
spin-polarized band-structure of Cu2Te20sBr2 and 
Cu2Te205Cl2 near the Fermi-level along various symme- 
try directions in the Brillouin zone. Calculations have 
been performed within the framework of the full poten- 
tial linearized augmented plane wave (LAPW) method 
based on the WIEN97 [6] code, the full potential min- 
imum basis local orbital code (FPLO) [7] and the lin- 
earized muffin tin orbital (LMTO) method based on the 
Stuttgart TBLMTO-47 code [8]. The band structures 
obtained by the three methods are in overall agreement 
with each other. 

The four narrow bands of width ^0.7 eV are well se- 
parated from the occupied low-lying valence band by a 
gap of ~ 0.25 eV and the high-lying unoccupied Te-p 
bands by a gap of ~ 2 eV. The bands near the Fermi- 
level do not contain essentially any Te-orbitals. They are 
of dominant Cu-3<i character (predominantly 3d x 2_ y 2 in 
the local frame of reference) as shown in the partial den- 
sity of states (DOS) in Fig. 3 with substantial admixture 
with Br(Cl) p, 01 and 02 p states as shown in Fig. 4. 

We can investigate these results further by analizing 
the electron-density. In Fig. 5 we present the electron 
density for Cu2Te20sBr 2 for one tetrahedral-unit. The 
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lobes of the Cu-3rf orbitals are oriented towards the near- 
est neighbor Br, 01, 01 and 02 ions, which form the dis- 
torted square surrounding of the Cu ion. The covalent 
bonding in between the Cu-3d orbital and the respective 
halogen-p and oxygen-p orbitals leads to the admixture 
of halogen and oxygen character at the Fermi-level, as 
seen in the respective DOS (Fig. 4). 

Comparing the band-structure of both compounds 
near the Fermi level (see Fig. 2) we note that the bands of 
both systems are quite similar and only differ in some de- 
tails that will translate in quantitative differences in the 
behavior of the effective model. The band dispersion, 
due to inter-tetrahedral matrix elements, is substantial 
(within the narrow bandwidth) along all three crystal- 
lographic directions indicating that the inter-tetrahcdral 
couplings are non- negligible in all three directions. 

The DFT-calculations yield, within the LDA or GGA- 
approximation, to four half-filled (metallic) bands. We 
performed a spin polarized calculation which leads to an 
antifcrromagnetic groundstate in agreement with the ex- 
periment, although the resulting band gap is too small. 
An additional local Coulomb-repulsion - not taken fully 
into account in the LSDA or GGA-approximation - will 
basically further shift these bands, enlarging the gap and 
resulting in the picture of lower- and upper-Hubbard 
bands in agreement with experiment. In the follow- 
ing, however, we did not do the simulation of the miss- 
ing Coulomb interaction in an LDA+C7 type calculation, 
which is in general useful but we do not expect from it 
new insights about the important couplings. 

Attempts to directly compute the exchange integrals 
by comparing the LSDA total energy differences for dif- 
ferent spin arrangements to that of a Heisenberg-like 
model faced two basic short-comings: (i) The resulting 
magnetic moments at the copper sites were quite dif- 
ferent for different spin arrangements. This leads to a 
strong bias of the inter-atomic exchange energies due to 
different intra-atomic contributions to the total energy, 
(ii) The important interactions in these compounds have 
a rather long-ranged nature as mentioned below. This 
would demand the calculation of large unit cells beyond 
the present computational capabilities. Therefore, we fo- 
cus this letter to accessing the corresponding transfer in- 
tegrals thereby pointing out the important interaction 
path-ways in these systems. 

Effective Model.- In order to quantify the results ob- 
tained from the ab-initio calculations in terms of hopping 
matrix elements, tij, we have employed LMTO-based 
downfolding [9] and tight-binding analysis on the band 
structure of these systems. The downfolding method con- 
sists in deriving a few-orbital effective hamiltonian from 
the full LDA or GGA hamiltonian by folding down the 
inactive orbitals in the tails of the active orbitals kept 
in the basis chosen to describe the low-energy physics 
of the system. This procedure naturally takes into ac- 
count the proper renormalization effect of integrated-out 



inactive orbitals in the effective interactions defined in 
the basis of the active orbitals. Using the real-space de- 
scription of the downfolded Hamiltonian one gets Hr = 

~ ( c ] c i + c \ c j) wnere Uj provides the effective 

hopping matrix elements between the active orbitals. In 
Table I we present the results for the most significant hop- 
ping matrix elements -shown schematically in Fig. 6- ob- 
tained from downfolding the full LMTO Hamiltonian to 
effective Cu-only Hamiltonian by integrating out every- 
thing except the d x 2_ y 2 orbital for each Cu atom in the 
unit cell to reproduce accurately the four narrow bands 
close to the Fermi energy. The first and second column of 
Table I show the hopping parameters for Cu2Tc205Br2 
and CU2TC2O5CI2 respectively. 

The predominant matrix elements consist of a set of 
nine different hopping parameters, some of them quite 
long-ranged and which couldn't be neglected in order 
to get a good description of the energy bands. Apart 
from the intra-tetrahedra t\ and ti parameters, the inter- 
tetrahedra t Xl t z , t ai t Cl t Vl td and t r are important (see 
Fig. 6). For instance, t c is needed in order to get disper- 
sion along the path T — M and t a and td are essential in 
order to get the correct behavior along the path X — A. 
The need of including such longer-ranged hoppings like 
t r is set by the renormalization process of the downfold- 
ing procedure [10]. We note, that long-ranged hopping 
matrix elements have proven to be essential for the de- 
scription of some related copper systems [11]. 

In order to investigate more in detail the nature of the 
Cu-Cu interaction paths and in particular the role of the 
halogen ion, we have also performed, a TB-downfolding 
analysis of the bandstructure of both systems by keep- 
ing the halogen X-p orbitals active in addition to four 
Cu-rf orbitals in the basis set (Cu+X downfolding, see 
third and fourth column of Table I). The results will be 
discussed in the following. 

TABLE I. TB-Downfolding parameters for the effec- 
tive Cu-Cu hoppings in meV units obtained by keep- 
ing active (i) the Cu-3d orbitals (first and second 
columns, Cu-downfolding) and (ii) both the Cu-3d with the 
halogen-3p/4p orbitals (third and fourth columns, Cu+X 
downfolding) . 
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Intra-tetrahedral couplings.- The nearest neighbor 
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(n.n.) intra-tetrahedral Cu-Cu transfer matrix element t\ 
is mediated by 01-ions located in between two Cu-ions 
and responsible for the Cu-01-Cu superexchange gen- 
crating the spin-spin coupling J\. This superexchange- 
path is evidenced in Fig. 5. The angle Cu-Ol-Cu 
is 107°/T09° for X = Br/CI. Assuming a similar 
crystal field, the larger superexchangc-anglc for X=C1 
results -following the Goodenough-Kanamori- Anderson 
rules [12]- in t\{X = CI) > ti(X = Br), in agreement 
with the downfolding results, see Table I. 

The intra-tetrahedral hopping matrix element ti cor- 
responds to the effective next nearest neighbor (n.n.n.) 
Cu-Cu overlapp. A comparison of the Cu- and the Cu+X 
downfolding result shows a substantial drop of t 2 for 
both systems when the halogen-orbitals are integrated 
out. This drop indicates that pathes of the type Cu-X- 
X-Cu are important for the effective ti matrix element 
and that they compensate to a certain extend the contri- 
bution of the direct Cu-Cu path (see Fig. 5) and of other 
possible pathes through oxygen ions. We note, that t 2 is 
vanishingly small in the Cl-compound. 

Inter-tetrahedral coupling.- The coupling t x via 
the 03 along the x— and y— direction in between two 
Cu4-tetrahedra (see Fig. 6 and Table I), is small, partly 
because the 03-weight is small at the Fermi-level. There 
is, however, a substantial inter-tetrahedral coupling t a 
and diagonal td within the .xy-planc mediated by the 
halogcn-p orbitals, (see Fig. 6 and Fig. 7). The role of 
the halogen ions for the inter-tetrahedral couplings can 
be quantified by analysing the Cu-downfolding and the 
Cu+X-downfolding results presented in Table I. Com- 
paring the two corresponding values (Cux and Cu+X) 
we observe that t a and td are nearly exclusively due to 
halogen-containing exchange pathes. 

The reason for this unusual large inter-tetrahedron 
coupling in Cu2Te205X2 is the large extension of the 
C\-3p and Br-4p wave functions which do not occur in 
cuprates containing only 0-2p orbitals (i.e. Cu04 units). 
In contrast to certain other cuprate compounds contain- 
ing halogen-ions like Sr2Cu02Cl2 [13], the situation for 
the Cu2Te205X2-compounds is unique. Here, the halo- 
gen is part of the covalent Cu-O-Cl(Br) network, whereas 
in the former mentioned compound family the halogens 
play only the role of an anionic charge reservoir [13]. 

A closer analysis of the electron density presented in 
Fig. 7 leads to two observations, which are in agreement 
with Cu+X downfolding results: (i) The Cl-3p orbital 
at the Fermi-level is more strongly distorted towards the 
Cu-3d orbital compared to the Br-4p orbital, indicating 
stronger copper-halogen coupling for X = CI than for 
X = Br (see arrows in Fig. 7) (ii) The Br-4p orbital is 
reoriented somewhat towards the n.n. Br-4p orbital indi- 
cating a stronger covalent Br-Br overlap. 

In addition to the inter-tetrahedral hopping processes 
along x and y there are substantial contributions to the 
inter-tetrahedral coupling terms t c , t v and t z along the 



^-direction, compare Fig. 6 and Table I. 

Discussion.- It has been recently argued [4] that 
Cu2Te20sBr 2 and Cu2Te205Cl2 differ only quantita- 
tively but not qualitatively in their magnetic properties, 
with Cu2Te205Br2 being somewhat closer to a quantum- 
critical phase transition. Here we present evidence from 
ab-initio calculations that the electronic properties of 
these two compounds are indeed close and that the inter- 
tetrahedral coupling is considerable. We note that a 
substantial coupling between Cu-tetrahedra is necessary 
in order to establish magnetic long-range order [4] . 

We have also presented a detailed analysis in terms of 
two different downfolding models (Cu and Cu+X) which 
reveal the important result that the halogen-ions are es- 
sential for the inter-tetrahedral exchange couplings. We 
find that four halogen orbitals, coupled by considerable 
(X-p)-(X-p) covalent bonding give rise to X 4 -rings which 
mediate long-ranged inter-tetrahedral couplings. These 
X4-rings are covalently coupled to the respective Cu 4 - 
tetrahedrons. 
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FIG. 2. Comparison of the band structure near the Fermi 
level between Cu2Te20sBr2 (full squares) and Cu2Te20sCl2 
(empty circles), with F = (0,0,0), X = (tt, 0, 0), 
M = (7r,7r,0), Z = (0,0, tt), R = (0,7r,7r) and A = (7r,7r,7r). 
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FIG. 4. Density of states for the (total) halogen-p and 
the Ol-p and the 02-p orbitals, for Cu2Te20sBr2 (left) and 
Cu2Te 2 5 Cl2 (right). 



FIG. 1. Crystal structure of Cu2Te205X2. The small green 
balls in the center of the bonds represent Cu, the medium 
sized yellow and large blue balls O and the halogen X, respec- 
tively. The atom not conected to a bond is 03. For simplicity 
the Te atom was droped. The inset presents an idealized 
Cu4-tetrahedron with the four corner-sharing Ol-Ol-02-X 
idealized squares and the Cu 2+ -ion in the center. 
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FIG. 5. Electron density for Cu2Te20 5 Br2 at the 0.2 
(e/A 3 ) isovalue for one tetrahedron. The dashed and da- 
shed-dotted lines indicate the four Br-Ol-Ol-02 distorted 
squares with Cu placed in the respective centers, they share 
Ol-corners. The solid line indicates the Cu-Ol-Cu exchange 
pathes (tj in Fig. 6). 




FIG. 6. Illustration of the hopping matrix elements t a 
(a = 1, 2, x, z, a, c, d) in between the Cu 2+ ions which are 
located in the center of the O3X distorted square (X=Br,Cl), 
indicated by the shaded regions (not to scale, in reality the 
distorted squares share Ol-corners). The filled circles denote 
the direction of the halogen-sites. The parameters t r and t v 
correspond to the hopping elements between Cu 2+ ions at A 
and B and at M and N respectively and equivalent positions. 




FIG. 7. Electron density for Cu2Te 2 5 Br2 (left) and 
Cu2Te 2 5 Cl2 (right) at the 0.05 (e/A 3 ) isovalue. for a projec- 
tion on the xy-plane in between four Cu-tetrahedra (compare 
with Fig. 6, left). The four p orbitals at the center are Br-3p 
(left) and Cl-4p (right). The respective arrows indicate the 
distortion of the X-p orbitals. 
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